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The  mean  or  expected  value  and  the  variance  of  the  absolute  value  of  correlation  between 
equal-length  independent  random  sequences  have  been  determined  and  compared  with  the  same 
measures  for  full-period  pseudonoise  (PN)  sequences,  which  have  applied  extensively  in  spread  - 
spectrum  radio  communications  systems.  The  magnitude  of  correlation  is  often  of  primary  interest 
in  such  applications,  since  a strong  negative  correlation  is  as  significant  as  a strong  positive  correla- 
tion; hence  correlation  detectors  exist  which  function  on  the  basis  of  the  magnitude  of  correlation. 
The  comparison  reveals  several  advantages  of  the  Gold  codes  which  are  due  to  their  inherent 
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mutual  orthogonality.  The  study  of  random  codes  permits  the  performance  evaluation  of  such 
detectors  during  periods  of  noise-only  input  and  also  provides  a long-needed  point  of  reference 
for  PN  system  design. 
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MEAN  AND  VARIANCE  OF  THE  CORRELATION 
MAGNITUDE  OF  RANDOM  AND  PSEUDONOISE  SEQUENCES 


INTRODUCTION 

The  mean  or  expected  value  and  the  variance  of  the  absolute  value  or  magnitude  of 
the  correlation  between  two  equal-length  independent  random  aequences  are  derived  for 
all  unique  bit-integral  relative  phase  conditions.  These  are  then  compared  with  the  same 
measures  for  the  full-period  correlation  between  pairs  of  certain  specific  types  of  pseu- 
donoise (PN)  sequences  which  have  been  applied  extensively  in  spread-spectrum  radio- 
communications systems.  The  study  of  independent  random  sequences  reported  here 
provides  a long-needed  point-of-reference  for  designers  of  spread-spectrum  systems,  be- 
cause the  autocorrelation  and  crosscorrelation  properties  of  candidate  PN  sequences  may 
now  be  compared  with  the  same  quantities  derived  for  fundamental  random  sequences. 
Correlation  detectors  exist  which  function  on  the  basis  of  the  nugnitude  of  correlation 
between  a receiver  input  and  an  internally  generated  reference  sequence.  The  treatment 
of  random  sequences  presented  here  permits  evaluation  of  the  performance  characteristics 
of  such  detectors  during  periods  of  noise-only  input. 

Design  studies,  validated  by  measurements  of  hardware  implementation,  have  proved 
that  differentially  coherent  phase-shift-keyed  (DPSK)  modulation  formats  permit  realiza- 
tion of  a relatively  high  communications  efficiency  for  digital  radio-communications  sys- 
tems when  interference  is  an  insignificant  problem.  An  important  motivation  for  the 
development  of  spread-spectrum  techniques  stems  from  the  need  to  mitigate  the  effects 
of  severe  electromagnetic  interference.  Spread -spectrum  systems  have  been  developed 
which  reduce  interference  significantly  and  approach  the  same  communications  effi- 
ciencies demonstrated  by  DPSK  systems. 

Spread-spectrum  systems  employ  an  encoding  digital  bit  stream  upon  which  base- 
bend  informational  data  are  impressed  as  a modulation.  The  baseband  information  can  be 
recovered  in  the  demodulation  process  only  by  correlation  with  a properly  synchronized 
replica  of  the  encoding  signal;  however  the  correlation  process  permits  a discrimination 
against  the  uncorrelated  intervening  interference.  Certain  reproducible,  deterministic, 
pseudorandom  sequences  which  provide  the  wideband  low-power-density  encoding  signals 
needed  for  such  systems  tend  to  timize  the  use  of  the  communications  capacity  offered 
by  a channel  with  noise  and  interference.  These  sequences  have  statistical  properties 
similar  to  noise  and  hence  are  known  as  pseudonoise  (PN)  sequences. 

The  usefulness  of  PN  sequences  as  codes  in  spread -spectrum  communications  depends 
on  both  their  individual  autocorrelation  and  their  joint  crosscorrelation  properties.  Auto- 
correlation properties  are  important  for  assuring  rapid  and  proper  synchronization  of  the 
incoming  signal  with  the  locally  generated  replica  of  the  PN  encoding  sequence.  To 
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perform  the  synchronization,  a scanning  or  time  search  between  the  incoming  signal  and 
the  local  PN  sequence  (often  generated  by  a binary  shift-register  code  generator)  is  re- 
quired. Sequences  with  off-peak  autocorrelation  values  comparable  to  the  main  correla- 
tion peak  are  undesirable,  because  they  tend  to  confuse  the  synchronization  process  and 
prevent  the  full  interference  discrimination  or  processing  gain  realizable  with  proper  syn- 
chronization. The  crosscorrelation  properties  of  the  sequences  are  of  considerable  signif- 
icance in  multiplexing  applications  in  which  many  communications  systems  must  operate 
in  a common  frequency  channel,  with  each  communication  link  employing  a different 
PN  sequence.  The  crosscorrelation  between  any  two  sequences  used  in  the  system  must 
be  sufficiently  low  to  provide  the  particular  amount  of  interference  isolation  required 
between  the  multiplexed  channels  for  satisfactory  operation. 

Prior  studies  by  others  have  considered  the  statistical  nature  of  only  the  correlation 
function  itself  for  various  sequences.  However  the  properties  of  the  absolute  value  of  the 
correlation  function  are  often  of  primary  interest  in  applications  such  as  spread-spectrum 
communications,  where  a strong  negative  correlation  is  equally  as  significant  as  a strong 
positive  correlation.  Gold  [1,  2]  has  developed  a method  of  selecting  certain  special 
maximal-length  linear  PN  sequences  which  he  used  in  pain  and  called  preferred  pairs.  A 
maximal-length  sequence  is  the  longest  which  can  be  repetitively  generated  by  a given  size 
code  generator.  The  least  upper  bound  of  the  absolute  values  (or  magnitudes)  for  the 
crosscorrelation  of  these  maximal-length  preferred  pain  is  better  (smaller)  than  for  any 
other  two  arbitrarily  chosen  maximals  of  the  same  length  (except  for  some  minor  excep- 
tions which  have  peaks  only  slightly  less  than  that  bound).  Because  these  special  codes 
were  originally  identified  by  Gold  and  are  most  often  used  in  sets  of  two,  they  are  com- 
monly referred  to  as  Gold  pain.  Further,  and  of  considerable  significance  for  practical 
applications.  Gold  [1,  2]  proved  that  these  same  low  crosscorrelation  magnitude  bounds 
exist  for  large  families  or  related  encoding  sequences,  known  as  Gold  code  families,  where 
each  member  of  a family  may  be  generated  by  forming  modulo-2  linear  combinations 
of  a Gold  pair  of  sequences.  Such  codes  have  the  same  period  as  the  Gold  pair  used  to 
generate  them.  The  bound  for  the  magnitude  of  crosscorrelation  between  any  two  mem- 
ben  of  a Gold  family  of  sequences  thus  can  be  characterized  by  just  the  bound  appro- 
priate for  the  Gold  pair  used  to  generate  the  family. 

Since  a statistical  analysis  of  the  properties  of  the  magnitude  of  correlation  for  all 
possible  phase  shifts  of  two  purely  random  sequences  apparently  has  not  been  developed 
and  presented  previously,  expressions  are  derived  for  the  expected  value  and  variance  of 
the  full-period  correlation  magnitude  values  for  each  relative  phase  relationship  between 
a synchronized  random  sequence  and  another  independent  code  of  the  same  length.  (If 
a binary  sequence  is  random  and  independent  of  another,  the  expected  correlation 
magnitude  between  the  two  is  independent  of  the  specific  nature  of  the  other,  because 
in  either  case  the  probability  of  matching  any  pair  of  bits  is  1/2.  Furthermore  the  ex- 
pected crosscorrelation  magnitudes  for  two  independent  random  sequences  is  therefore 
the  same  as  their  individual  expected  autocorrelation  magnitudes,  for  the  same  reasons.) 
Random  code  lengths  corresponding  to  maximal-length  PN  sequences  have  been  selected 
for  study  to  facilitate  a direct  comparison  of  their  respective  characteristics.  The  analysis 
provided  by  this  study  is  applicable  to  bit-synchronized  data  with  negligibly  small 
Doppler-frequency  effects.  The  crossspectral  correlation  characteristics  which  are  discussed 
in  Figs.  3 and  4 for  nonsynch ronized  data  with  significant  relative  Doppler-frequency 
shifts  are  not  treated  here. 
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DERIVATION  OF  THE  EXPECTED  VALUE  FOR  THE  MAGNITUDE  OF 
THE  CORRELATION  BETWEEN  TWO  RANDOM  SEQUENCES 

For  time-dependent  periodic  functions  f(t)  and  g{t)  of  period  T,  the  normalized 
crosscorrelation  R(r)  between  them  is 


R(t) 


f(t)git  + T)dt. 


(1) 


The  functions  f(t)  and  g(t)  of  the  continuous  variable  t are  now  replaced  with  the  func- 
tions a(i)  and  b(i)  of  a discrete  variable  i (which  is  a common  form  for  describing  se- 
quences), so  that  the  crosscorrelation  may  be  expressed  in  the  analogous  discrete  form 


1 L 

m - l£o(i)b(t+j). 


(2) 


I 

< 


where  the  a(i)  and  b(i  + ;')  values  are  taken  from  the  set  {+ 1,  -1}  and  L is  the  number 
of  discrete  elements  in  the  period.  Hence 

a(«)b(«  +;)  = +1,  if  a(i)  - b(i+j), 

- -1,  if  <i(i)  * b(i+j). 


(3) 


If  n represents  the  number  of  agreements,  then  L - n will  represent  the  number  of 
disagreements;  and  R(J)  can  be  written  as 


R(J)  m L 


2]  o(/)b«  *j)  + ^ 

i i 

L •(/) -*(<♦»  o U)+b(i+j) 


o(l)b(i  *j) 


i 

<K0"  8(1  ♦/) 


L 

i 


(4) 


(5) 


(6) 


The  expression  in  brackets  in  (6)  is  the  form  given  by  Gold  (1,  2]  for  the  correlation, 
namely,  the  number  of  agreements  minus  the  number  of  disagreements.  The  multiplicative 
factor  (1/L)  normalizes  the  correlation.  Following  Gold  [1,  2]  0 will  be  used  to  represent 
the  normalized  discrete  correlation  9(n),  (where  the  discrete  variable  n here  represents  the 
number  of  agreements  between  the  two  sequences,  rather  than  a particular  phase  shift  of 
the  second  sequence).  Thus 


0(n)  - (2 n/L)  - 1. 


(7) 
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Given  two  sequences  of  the  same  length,  if  at  least  one  sequence  is  random  and 
independent  of  the  other,  the  probability  of  an  agreement  between  any  two  correspond- 
ing bits  is  1/2.  Hence  the  probability  of  n agreements  arranged  in  a particular  order  is 
(l/2)n.  Similarly  the  probability  of  L - n disagreements  occurring  in  a particular  order  is 
(l/2Vt _n.  With  use  of  the  binomial  coefficient,  namely, 

(L)‘  —kL~  - - L!  -(  L ) (8) 
W n!(L-n)!  (L-n)'.n'.  \L-n)’  w 

which  is  the  number  of  different  ways  of  obtaining  n agreements  and  L - n disagreements, 
the  probability  of  getting  exactly  n agreements  in  any  order  is 


(W 

m- 


The  expected  value  E(  |0|),  or  mean  n,  of  the  absolute  value  of  the  correlation  is  the 
sum  of  each  of  the  possible  values  of  |0(n)|,  weighted  by  the  corresponding  probability 
of  occurrence  for  each  case.  Thus 


u * £101  = 2]  • 

n-0 


5't  ; 


■ & £*  ■ <)■ 


By  definition  of  the  absolute  value. 


II  - 1 - — , ifO<  n < K, 


m - - 1.  if  JT  + 1 < n < L, 


where 


Lodd, 


L even. 
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Thus  (10)  can  be  written 
'1\L 


Let  m m L - n;  then,  by  use  of  (8),  the  second  summation  in  (13)  becomes 

■t ('”)!-)' 


(13) 


0 

I 

which,  after  the  index  m is  replaced  with  n,  becomes 


1, (?•')(.')■  5("?)0 


-£(-30  ■£('•?» 


L even, 


L('-m 

n-0 


or  even. 


(14) 


The  summand  for  n - JT  ♦ 1 is  neglected  for  L even  because  this  summand  is  0;  that  is,  if 


l 


n - K + 1 


then 


(18) 


i 
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Substitution  of  (14)  in  (13)  gives 

"ffl‘  [#•?)(!)] 

■(T©:)- £?«)]■ 

(16) 

where  in  the  second  summation  the  lower  index  becomes  1 instead  of  0 because  n is  a 
multiplicative  factor. 

Hie  binomial  expansion 

(a+by  = 

becomes,  for  a = b m 1/2, 

(17) 

■ • goer 

•ergo- 

(18) 

i 

which  can  be  written 

* ■ £0  • 
i-0'  ' 

(19) 

In  the  notation  desired  here,  (19)  is  written  as 

21  ■ £(») 
/i-0 

v a 

-£$♦  tft- 

♦ 

(20)  | 

which,  by  use  of  (12)  and  m - L - n,  can  be  written  as 
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11 

< 

( 

i 
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m-tn  - 


(21) 


Then,  if  the  arbitrary  index  m is  replaced  by  n and  the  order  of  summation,  is  reversed, 
(21)  becomes 


or 


K-(L- 1)/2 

Z 

n- 0 

K-(X.-2)/2 

Z 

n-0 

Also,  from  the  definition  of  the  binomial  coefficient. 


Km(L~\  )/2,  . . 

*'■-*  Z („)'t,,dd- 

n»0  ' ' 

_ oKm{^r)l2(L\  ( L \ t 

£ U/  \L/2/,L 


even. 


(L-l)! 


(n-l)![(L-l)  - (n  - 1)] ! 


Use  of  (23)  and  (24)  in  (16)  yields 


(22) 


(23) 


(24) 
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■ri-'-#::)]-- 

•(ri’-'-St)'#:.')]- 


is  introduced,  then 


where 


M = n - 1 


L-  3 

K - 1 = — , L odd, 


= — - — , L even. 

When  L is  even,  L - 1 is  odd,  and  vice  versa;  therefore  if  (23),  is  applied,  to  the  case  in 
which  L is  replaced  by  L - 1,  it  can  be  written  as 


(L-3)/2/r 

(Y)- 


L odd. 


(L- 2)12 


SI*;1). 


L even. 
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For  L even,  the  expression  in  brackets  may  be  simplified  using  (8): 
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In  expanded  form  p ■ 1 for  L m 1 and  for  L > 2, 


(34) 


even. 


By  use  of  Stirling’s  approximation  for  factorials  namely, 

m!  e'mrnm(2irm)1/2, 


(35) 


the  following  approximate  evaluations  for  (33),  valid  for  L > > 1,  can  be  obtained; 

w (1/2 )w  [e-wmw(2PM)1/23  % (lV12 

Mrwrr  ' 

where 

2 < m “ L - 1,  L odd, 

• L,  L even.  (36) 

This  approximate  relationship  is  particularly  useful  when  L is  very  large,  because  the 
direct  evaluation  of  (34)  takes  excessive  computer  time.  Table  1 presents  a test  program 
for  comparing  (33),  or  (34),  with  (36)  for  several  sequences  of  period  L * 2N  - 1,  where 
N is  a positive  integer  corresponding  to  the  number  of  stages  in  a binary  shift-register 
code  generator.  The  results,  shown  in  Table  2,  indicate  that  (36)  becomes  a progressively 
better  approximation  as  N (and  L)  increases  and  that  the  error  is  negligible  (less  than 
0.0003  dB)  for  L - 10,000  or  more. 


DERIVATION  OF  THE  VARIANCE  OF  THE  MAGNITUDE  OF  THE 
CORRELATION  BETWEEN  TWO  RANDOM  SEQUENCES 

The  variance  o*  of  the  absolute  value  of  the  correlation  about  the  mean  value  is  de- 
fined as  the  sum  of  the  squares  of  the  deviations  from  the  mean,  weighted  by  the  prob- 
ability of  occurrence  for  each  case.  Ulus,  by  use  of  (7)  and  (10)  and,  by  definition, 

L 

X]  P(«)  “ 1. 

n-0 
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Table  1 — FORTRAN  Program  to  Validate  the  Approximation  Given  by  (36) 


OOiOO  PR06RAH  STIRL <0UTPUT> 

001 10  PRINT'*  N L 

00120*-'*  ♦ 

00130  DO  10  N-3»17 
00140  L-2**N-1 
00130  CALL  AVT  <L»E> 

00160  CALL  XAVT  <L*  T> 

00170  DDB-20.  ♦ALD61 0 CTVE) 

00180  PRINT  20»N»L»E'T'DDB 
00190  10  CONTINUE 
00200  20  PORMAT <2I8»2F10.6*E12.3> 
00210  END 

0022  0C 

0023  0C 

00240  SUBROUTINE  AVT<L»E> 

00250  J-L 
00260  E-L^2 
00270  K-E*2. 

00280  IF<K.EQ.L> J«J*1 
00290  K-0 
00300  J«  <J-1) '2 
00310  E-l. 

00320  00  10  I-1»J 

00330  X-K-K-U  S Y-K-K+l 

00340  E»E*X^Y 

00330  10  CONTINUE 

00360  RETURN 

00370  END 

00380C 

0039 OC 

00400  SUBROUTINE  XAVT<L»T> 

00410  B* 0. 6366 1 9772365733 
00420  J-L 
00430  A-L/-2 
00440  K-A*2. 

00430  IF <K.EQ.L> J*J*1 
00460  A-J-l 
00470  T-SQRT  <B^A) 

00480  RETURN 
00490  END 


TRUE 

VALUE 


APPROX. 

VALUE 


ERROR*' 
IN  DB*'/' 
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Table  2 — Computed  Result*  of  the  Validation 
Program  of  Table  1 


N 

L 

TRUE 

APPROX. 

ERROR 

value 

VALUE 

IN  D» 

3 

7 

.312500 

.323735 

. 36029E+00 

4 

15 

.209473 

.213244 

. 13497E+00 

3 

31 

.144464 

.145673 

. 72369E-01 

6 

63 

. 1 00924 

.101331 

. 33022E-01 

7 

127 

.070940 

. 071081 

. 17234E-01 

8 

255 

.050014 

. 050064 

. 83491E-02 

9 

311 

.033314 

.035331 

. 42578E-02 

10 

1023 

. 024952 

. 024938 

. 21247E-02 

11 

2047 

.017637 

.017640 

. 10613E-02 

12 

4095 

. 012469 

.012470 

. 53 04 0E- 03 

13 

8191 

. 008816 

.008817 

. 26314E-03 

14 

16383 

. 006234 

. 096234 

. 13255E-03 

IS 

32767 

. 004406 

. 004408 

. 66273E-04 

16 

63533 

.003117 

.003117 

. 33136E-04 

1? 

131071 

. 002204 

. 002204 

. 16570E-04 

the  variance  is 

L 

a2*  J^mn)\-H)2p(n) 

n-0 

L L L 

- 2]  1®(«)12P(»)  - 2m  W")IP(«)  ♦ P2J^P(n) 

»■ 0 n*0  n-0 

L 

- £«*<n)p(n)  - 2n(m)  + M2(l) 

L 

- £«2(n)P(n)  - m2 

«-o 


(37) 


(38) 
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where  the  index  has  been  shifted  from  n*0ton*lto  account  for  the  zero-value 
summand.  By  use  of  (8)  the  last  summation  becomes 


y*  n / L\  _yi  nLl 

hL\nihLn[(L-n)- 


From  (19)  and  (8) 


r - 1)' 

L>x  (n  - 1)![£  - 1 - (n  - 1)J ! 

£(1:3 

t-Z-tfi'-l) 

‘*£«) 


2*  - y -L! 

n!(L-n)! 


-**£o 


i 
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When  applied  to  the  case  L ■ P - 1 and  q * n + 1,  (41)  becomes 

1 + 


-‘-Sr;1) 

which,  by  an  exchange  of  arbitrary  parameters,  can  be  rewritten  as 

2l'x  - 1. 


£(;:;) 


Substitution  of  (43)  into  (40)  gives 


in.)-*-1- 

n*  1 ' 


which,  when  substituted  into  (39),  gives 

\i-2i 


1 


where  J • L - 1 and  A#  ■ n - 1.  Then 


-•sir 


- 1 - M2. 


From  (44)  and  (19),  (45)  simplifies  to 
l/l^*1 


08  " l(l)  1 + " 1 _ **» 

which  further  simplifies  to 


(42) 


(43) 


(44) 


(45) 
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or  finally 


j2  = 


(J+2)  - 1 - /i2 
I 1 ~ "2 


L 

L + 1 


(46) 


DETERMINATION  OF  THE  EXPECTED  VALUE  FOR  THE  MAGNITUDE  OF 
CORRELATION  BETWEEN  TWO  GOLD-PAIR  SEQUENCES 

To  obtain,  for  comparison,  the  distinguishing  characteristics  between  purely  random 
and  PN  Gold  codes,  the  expected  value  and  variance  of  the  magnitude  of  crosscorrelation 
for  a preferred  pair  of  PN  Gold  sequences,  over  the  various  possible  code-phase  combina- 
tions, is  also  determined.  Again,  as  in  (10),  where  now  L = 2N  - 1,  the  mean  or  ex- 
pected value  is 


■ L l*M-  • 

i 


But  now 


(47) 


Pi  = fi/(2*  - 1) , (48) 

where  /j-  is  the  frequency  or  number  of  occurrences  of  a particular  value  of  |0(|  over  the 
period.  Since  0,  is  a deterministic  function  for  Gold  pairs,  the  fs  can  be  uniquely  determined 

Gold  (5J  has  shown  that  the  correlation  for  a Gold  pair  when  N is  odd  can  only  be 
one  of  three  specific  values,  namely. 


[-2<N*i)/2  - 11/(2* -1), 

(49a) 

[2<N*l)/2  - l]/(2*-l), 

(49b) 

-1/(2* -1), 

(49c) 

with  frequencies  of  occurrence 
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fx  = 2""*  - 2<"'3)/2  , (50a; 

f2  = 2"‘2  + 2(^-3)/2  , (50b; 

f3  = 2"’1  - 1,  N odd.  (50c; 

Substitution  of  (48),  (49),  and  (50)  into  (47)  gives  for  N odd 

= [2<"+1>/2  + Jj  j 2n -2  _ 2<tf-3)/2] 

# (2"  - l)2 

^[2(**1)/2  + ij  j2N.2  + 2(Jv.3)/2j  + aw-i  - 1 
(2"  - l)2  (2"  - l)2 

p,  = [2<3"_1>/2  - 2<"_1>/2  ♦ 2"’l  - l]/(2"  - l)2,  N odd.  (51 

Sherman  [3]  provides  the  distribution  for  N even,  where  N = 6 + 4X  and  X = 0, 

1,  2 (Pairs  of  maximal  length  sequences  with  three-level  crosscorrelation  values  do 

not  exist  for  N = 0,  4,  8,  12 ) In  this  case  the  three  correlation  values  and  frequen- 

cies of  occurrence  are 

6X  = [-2<"*2)/2  _ j]/(2"-l),  (52a 

02  = [2<w*2)/2  - l]/(2"-l),  (52b; 

03  = -1/(2"  - 1) , (52c] 


fx  = 2" -3  - 2<"-<)/2  , 
f2  = 2" -3  + 2<"-4>'2, 

f3  = 3(2"-2)  - 1 , N = 4X  ; X = 0,  1,  2,  3, 
Substitution  of  (48),  (52),  and  (53)  into  (47)  gives 


• *2"  - n* 


(2"  - 1)2 

♦ t2*'a  + [2(N.„/a _ n + 

(2" -l)2  1 ' (2" -1)2 


H,  . [2<3"-*)/2  - 2("-*)/2  ♦ 3(2"-*)-  11/(2" -!)2,  N = 6 ♦ 4X;  X-  0,1,2 
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DETERMINATION  OF  THE  VARIANCE  OF  THE  MAGNITUDE  OF  THE 
CORRELATION  BETWEEN  TWO  GOLD-PAIR  SEQUENCES 

The  variance  of  the  magnitude  of  the  correlation  function  for  PN  Gold  pairs  can  be 
obtained  in  a similar  manner.  Thus,  for  N odd,  by  use  of  (38),  (48),  (49),  and  (50), 


°2  = £ Oft  - M,2 


[2(N+l)/2  + 


l2*-2  - 2<*-3>/2l 


(2*  - l)3 

+ 7 (2^-2  + 2<jv‘3)/2J  + ~ 1 _ M2 

(2W  - 1 )3  (*>N  - 1 »3  3 ’ 


(2*  - l)3 


22N  _ 2N  - 1 

(2N  - l)3 


- m;  . N odd  . 


For  N = 6 + 4X,  A = 0,  1,  2,  ....  by  use  of  (38),  (48),  (52),  and  (53), 
2 [2!iN*2)l2  + u2 


(2*  - l)3 


[2 N-8  - 2{N-4)I2  J 


(55) 


f2(N+2)/2  _ jl2  3(2n~2)  - 1 

+ — J [2*'3  + 2<w-0/2]  + L_L_  u2 


(2^-  l)3 


(2*-l)3 


(56) 


The  result  is  that  (55)  and  (56)  are  of  the  same  form. 


THE  CROSSCORRELATION  AND  OFF-PEAK  AUTOCORRELATION 
UPPER  BOUND  FOR  THE  MEMBERS  OF  GOLD  CODE  FAMILIES 

As  previously  stated.  Gold  [1,2]  proved  that  the  maximum  magnitude  of  the  cross- 
correlation  which  can  occur  for  a preferred  pair  of  maximal-length  sequences  is  also  the 
maximum  crosscorrelation  magnitude  or  bound  that  will  be  achieved  for  any  combination 
of  two  members  from  the  family  of  sequences  known  as  Gold-family  codes.  A method 
of  generating  a Gold  family  of  sequences  is  diagrammed  in  Fig.  1.  Code  A and  code  B 
are  maximal-length  linear  sequences  which  constitute  a preferred  or  Gold  pair.  Each  has 
a period  of  length  2N  - 1,  for  N- stage  digital  shift  registers.  A member  C,  of  the  Gold 
family  corresponding  to  codes  A and  B can  be  formed  by  the  modulo-2  addition  of  code  A 
with  a particular  phase  of  code  B.  The  sequences  formed  in  this  manner  belong  to  the 
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Fig.  1— A method  of  generating  member*  of  ■ Gold  family  of  code* 
auch  that  each  member  ha*  a bounded  croeacorrelatioa  magnitude  be- 
tween itaelf  and  any  other  member  within  the  family.  r(  is  an  integral 
number  i of  clock  cycles  of  time  delay  for  each  member  C,  of  the 
family,  where  i ■ i,  2,  ....  2N  - 1.  Code*  A and  B are  selected  (Gold- 
Pair)  N-stage  shift-register-generated  maximal-len7*h  linear  sequences 
which  arc  also  member*  of  the  Gold  code  family. 


nonmaximal  class  of  codes;  however,  the  two  maximal-length  sequences  used  to  generate 
the  other  members  of  the  family  are  themselves  included  in  the  family.  Since  there  are 
2n  - 1 possible  phases  of  code  B,  there  are  2*  * 1 members  in  the  Gold  family. 

Golomb  [6]  gives  the  off-peak  normalised  autocorrelation  values  for  a maximal-length 
linear  sequence,  such  as  used  in  generating  a Gold  code  family  (such  as  code  A or  code  B), 
which  are  the  same  low  value  (the  negative  reciprocal  of  the  sequence  length),  irrespective 
of  the  relative  code  phase.  Thus  the  off-peak  normalized  autocorrelation  magnitudes  for 
either  code  A or  code  B are 

|0|  - 1/(2* -1).  (57) 

The  off-peak  normalized  autocorrelation  magnitudes  for  the  other  members  of  the  Gold 
family  are  not  as  ideal;  however  Gold  has  stated  in  a private  communication  that  the 
off-peak  autocorrelation  values  for  all  codes  in  a Gold  family  are  bounded  in  magnitude 
and  that  this  upper  bound  for  any  one  family  member  is  the  same  as  the  crosscorrelation 
amplitude  bound  for  any  two  family  members.  This  bound  is  known  as  the  Gold  bound. 
The  largest  crosscorrelation  magnitude  of  the  Gold  pair  used  in  generating  the  family  can 
thus  be  used  to  characterise  the  crosscorrelation  off-peak  and  autocorrelation  bounds  for 
the  whole  family  of  sequences. 

The  Gold  bound  for  the  normalised  croascorrelation  and  off-peak  autocorrelation 
ms#iiti»rVe  over  the  full  period  for  a Gold  pair  (or  for  any  two  members  of  the  associ- 
ated family  of  codes)  is  given  by 

|0|  < l#*4*>'*  ♦ lJ/(2*  - 1) , Ilf  odd  , 

< [2 (*♦*>/*  ♦ 11/(2*  - 1) , N even  . (58) 

With  N odd,  this  upper  bound  is  a sharp  bound  in  the  sense  that  some  relative  phase  shift 
of  the  pairs  of  sequences  can  be  found  that  provides  a magnitude  of  correlation  which  is 
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equal  to  the  bound.  With  N even,  however,  there  are  some  non-Gold-pair  sequence  com- 
binations which  yield  maximum  values  for  |0|  which  are  just  slightly  less  than  the  bound 
indicated  in  (58). 


A COMPARISON  OF  PSEUDONOISE  GOLD-PAIR  CODES  AND  INDEPENDENT 
RANDOM  SEQUENCES 

A FORTRAN  computer  program  has  been  devised  to  evaluate  the  several  expressions 
providing  the  mean,  variance,  and  bounds,  where  applicable,  for  two  random  sequences 
and  Gold-pair  sequences  of  lengths  L = 2N  - 1.  The  program  listing  is  provided  as 
Table  3.  Table  4,  Fig.  2,  and  Fig.  3 give  the  results  as  computed  by  the  program. 

The  program  is  based  on  equations  (34),  (36),  (46),  (51),  (54),  (55),  (56),  (57),  and  (58). 
The  approximate  relation  given  in  (36),  based  on  Stilling’s  approximation,  is  used  when 
the  length  is  greater  than  10,000  in  order  to  avoid  excessive  computation  time. 

The  term  Gold  bound  as  used  here  is  synonymous  with  bounds  on  the  magnitudes  of 
the  normalized  crosscorrelation  of  a Gold  pair  (or  for  the  autocorrelation  or  crosscorrelation 

N 

3 4 S < 7 8 9 10  11  12  13  14  16  IS  17  18  19  20  21  22  23 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


Fig.  2— Tha  bounds  and  expected  magnitudes  of  tha  normalised  full-pa riod  autocorrelation 
and  cross conslation  of  Ooid  family  cods  sequences.  (0  dB  represents  the  largest  possible 
correlation  nra^iitude  which  occurs  when  there  is  either  total  agreement  or  disagreement 
among  corresponding  bits  of  the  sequences.) 
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Table  3 — FORTRAN  Program  for  Generating  Table  4 


00100  PROGRAM  BCORR  (OUTPUT) 

00110  PRINT  200 
00120  SQ2-SQRT(2.) 

00130  DO  10  N-3.23 
00140  o«N 
00130  H-N-/2 
00160  H*H*2 . 

00170  F-N/4 
00180  F-F*4. 

00190  L*-l *2**N 

00200  P-L 

00210  £-0.5*(l.*N> 

00220  T*2. ♦♦£ 

00230  IF(6.EQ.H)T«T*SQ2 
00240  T-T*l. 

00230  GB-TxP 

00260  6BD*20.*AL0610(GB> 

00270  CALL  AVTHETA  <L * RM) 

00280  RMD»20. ♦ALDG1 0 (RM) 

00290  RV-1"L-RM*RM 

00300  U2SV-20. ♦AL0G1 0 (RM*2. ♦SORT (RV) ) 

00310  IF( G.£Q.H)CALL  EVEN(N.GM.GV) 

00320  IF (G.NE.H)CALL  0DD(N.6M.GV) 

00330  GMD»2 0 . *AL06 1 0 (GM> 

00340  U1SV-20. ♦ALOGl 0 (GM+SQRT (GV)  ) 

00330  6AD»20.*AL0610(1./'P> 

00360  PRIMT  210th 
00370  IF (F.EQ.G) PRIMT  230 
00380  IF (F.NE.6)PRINT  240 

00390  IF (F. ME. G) PRINT  250.P.6B.GM.GV.RN.RV.6AD.6BD.6MD.U1SV.RMD.U2SV 
00400  IF (F.EQ.6)PRINT  233.P. GB. RM. RV, GAD. GBD. PHD. U2SV 
00410  10  CONTINUE 
00420  PRINT  260 

00430  200  FORMAT (18X.*EXPECTED  (AUTO^CROSS)  CORRELATION  MAGNITUDES*. ✓ . 

0 044 0*  1 5X » *BE  TWEEN  PAIRS  OF  MAXIMAL  GOLD  AMD  RANDOM  SEQUENCES*."/'. 
00430*2X. *N*. 5X» *C0DE  ♦. . 3 (7X. *GOLD*> . 2 (3X. *RANDOM*> . " 

0046 0*6X. *LEN6TH*. 7X.*B0UND*. 7X. *MEAN*. 3X» ♦ VARIANCE*. 7X. 

00470**MEAN*. 3X. ^VARIANCE*." . 3X. *60LD  AUTO-DB*. 7X. 

00480**DB*. 9X»  *DB*»5X» *U1 SV (DB) *. 7X. 

00490**DB*.  3X»  *U2SV  (DB)  *.  2H**"> 

00300  210  FORMAT (IX.  12. '*) 

00310  230  FORMAT (1H*.~> 

00520  240  FORMAT (♦  ♦ "> 

00330  250  FORMAT  (F8.  »3X»F9.  7.4F1 1. 7".  FI  7. 2.F9. 2.F  1 1 . 2.F  13. 2«  F9. 2.  F 1 3. 2. /") 
00340  253  FORMAT (F8. . 3X. F9.  7. 22X. 2F 1 1 . 7". F 1 7. 2.F9. 2. 22X.F1 1 . 2, F 1 3. 2") 
00330  260  FORMAT (". IX. IN*.*  IF  N IS  A MULTIPLE  OF  4.*. 

00360**  THE  NUMBER  OF  AVAILABLE  GOLD  PAIRS  IS  LIMITED.*. "* 

00370*1X,2H**.*  U1SV  ■ UPPER  1-SIGMA  VALUE*. 

00380**!  U2SV  - UPPER  2 -SIGMA  VALUE**". 

00390**  NOTE«  CORRELATION  BOUND  • ♦. 32M(2**( (N*!)/?) *1) ✓ (2**N-1)  . N ODD". 
00600*23X.33H>  (2**  ( (N*2)  '2> ♦ 1 ) ✓ (2**M- 1 ) . N EVEN.". 

0061 0*1 8X.  *LEN6TH  - 2*.2M**.*N-1*""> 

00620  END 
00630C 
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Table  3 (Continued)  — FORTRAN  for  Generating  Table  4 


0064 0C 

00630  SUBROUTINE  RVTHETfl <L . RH> 

00660  J-L  * RH-L/2  * K»RH*2. 

00670  IF<K.EQ.L>J-J+1 

00680  IF<L.LT. 10000>GO  TO  5 

00690  RH-SORT  <0. 6366 1 97723657493^  < J- 1 > > 

00700  GO  TO  20 

00710  3 CONTINUE 

00720  K-0 

00730  J-<J-l>/2 

00740  RH-1. 

00730  DO  10  I-l.J 

00760  X-K-K+l  S Y-K-K+l 

00770  RH*RH*X/Y 

00780  10  CONTINUE 

00790  20  CONTINUE 

00800  RETURN 

00810  END 

00820C 

00830C 

00840  SUBROUTINE  ODD<N*GM*GV> 

00830  FI *2. •♦<-N> 

00860  F2-1.-F1 

00870  F3-2. ♦♦ <-  <1*N> ^2> 

0088  0 F4-F3-2.  ♦♦  <-  <3»t1+ 1 > ^2)  +2.  ♦♦  <- 1 -N) 
00890'*— 2. <-2*N> 

00900  GH*  <F4^F2) /-F2 
00910  F5-2 . <N-1 ) '2> 

00920  F6-2. <- <N»3> ^2> 

00930  GV*  <F3+F 1 > **2* < 0 . 23-F6> 

00940  6V-GV+  <F3-F1> •♦2*<0. 23*F6> 

00950  6V-6V+2. ♦♦ <-2*N- 1 ) -2. <-3*N> 

00960  GV- < <GV/F2> /F2> /F2-6H*6N 

00970  RETURN 

00980  END 

0099 0C 

01000C 

01010  SUBROUTINE  EVEN<N» GK» GV> 

01020  Fl-2. 

01030  F2-1.-F1 

01040  F3-2. ♦♦ <- <2*N> ^2) 

01030  F4-F3-2. <- <3*N«-2> '2>  +3. *2.  ♦♦  <-2-N> 
01 060*-2. <-2*N> 

01070  GH« <F4xF2> ^F2 
01080  F3-2.**<-<N-2>/'2) 

01090  F6— 2. <H*4> '2) 

01100  GV*  <F5>F 1 ) **2* <0.1 25-F6> 

OHIO  GV*6V-MF3-F1>**2*<0.  125+F6) 

01120  GV-6V+3. *2. ♦♦ <-2*N-2> -2. ♦♦ <-3*N> 
01130  GV- < <GV/'F2>  ^F2>  ^F2-6N*GM 
01140  RETURN 
01130  END 
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Table  4 — Expected  (Auto/Cross)  Correlation  Magnitudes  Between 
Pairs  of  Maximal  Gold  and  Random  Sequences 


N 

CODE  GOLD 

LENGTH  BOUND 

GOLD  AUTO-DB  DB 

GOLD 

MEAN 

DB 

GOLD 

VARIANCE 

U1SV(DB> 

RANDOM 

MEAN 

DB 

RANDOM 

VARIANCE 

U2SV  <DB> ♦♦ 

3 

7 .7142837 

. 3469388 

. 0399833 

3125000 

. 0452009 

-16.90  -2.92 

-9.  19 

-5.24 

-10. 10 

-2.64 

4^ 

15  .6000000 

• 

2094727 

. 0227879 

-23.52  -4.44 

-13.58 

-5.23 

5 

31  .2903226 

. 1446410 

.0123441 

1444644 

.0113861 

-29.83  -10.74 

-16.79 

-11.84 

-16.80 

-8.92 

6 

63  .2698413 

. 0753338 

.0104458 

1009237 

. 0056874 

-35.99  -11.38 

-22.46 

-15. 01 

-19.92 

-11.98 

7 

127  .1338383 

. 0668981 

.0034602 

0709403 

. 0028415 

-42.08  -17.47 

-23.49 

-18. 01 

-22.98 

0 
in 

^4 

1 

8* 

235  .1294118 

. 

0500145 

. 0014201 

-48.13  -17.76 

-26. 02 

-18. 04 

9 

511  .0645793 

. 0322877 

.0009183 

0353136 

. 0007099 

-54.17  -23.80 

-29.82 

-24. 07 

-29. 04 

-21. 05 

10 

1 023  . 0635386 

. 0163732 

.0007104 

0249522 

. 0003549 

-60.20  -23.94 

-35. 72 

-27.33 

-32. 06 

-24. 06 

11 

2047  .0317538 

. 0158768 

.0002367 

0176374 

. 0001774 

-66.22  -29.96 

-35.98 

-30. 10 

-35. 07 

-27. 08 

12* 

4095  .0315018 

• 

0124692 

. 0000887 

-72.25  -30.03 

-38. 08 

-30. 09 

13 

8191  .0157490 

. 0078745 

.0000601 

0088163 

. 0000444 

-78.27  -36.05 

-42. 08 

-36. 12 

-41.09 

-33.10 

14 

16383  .0156870 

. 0039523 

. 0000454 

0062339 

. 0000222 

-84.29  -36.09 

-48. 06 

-39.42 

-44.  10 

-36. 1 1 

13 

32767  . 0078433 

. 0039216 

.0000151 

0044079 

. 0000111 

-90.31  -42.11 

-48. 13 

-42. 14 

-47. 12 

-39. 12 

16* 

65535  . 0078279 

0031 168 

. 0000055 

-96.33  -42.13 

-50. 13 

-42. 13 

17 

131071  .0039139 

. 0019570 

.0000038 

0022039 

. 0000028 

-102.35  -48.15 

-54.  17 

-48. 16 

-53. 14 

-45. 14 

18 

262143  .0039101 

. 0009794 

.0000029 

0015584 

. 0000014 

-108.37  -48.16 

-60. 18 

-51.47 

-56. 15 

-48. 15 

19 

524287  .0019550 

. 0009775 

.0000010 

0011019 

. 0000007 

-114.39  -54.18 

-60.20 

-54. 19 

-59. 16 

-51.16 

20* 

1048575  .0019541 

• 

0007792 

. 0000003 

-120.41  -54.18 

-62. 17 

-54. 17 

21 

2097151  .0009770 

. 0004885 

.0000002 

0005510 

. 0000002 

-126.43  -60.20 

-66.22 

-60.21 

-65.  18 

-57. 18 

22 

4194303  .0009768 

. 0002443 

.0000002 

0003896 

. 0000001 

-132.45  -60.20 

-72.24 

-63.52 

-68. 19 

-60. 19 

23 

8388607  . 0004884 

. 0002442 

.0000001 

0002755 

. 0000000 

-138.47  -66.22 

-72.25 

-66.23 

-71.20 

-63.20 

♦ IF  N IS  0 MULTIPLE  OF  4.  THE  NUMBER  DF  AVAILABLE  GOLD  PAIRS  IS  LIMITED. 


♦♦  UISV  - UPPER  1-SIGMA  VALUE!  U2SV  « UPPER  2-SIGMA  VALUE 

NOTE!  CORRELATION  BOUND  - <2**< <N*1> ✓£> *i> / <2**N-1> . N ODD 

- (£**  < <N*2>  ^2>  ♦ 1 > •''  <2**N- 1 > t N EVEN 

LENGTH  ■ 2**N-1 
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N 


3 4 S 6 7 8 9 10  11  12  13  14  15  16  17  18  19  20  21  22  23 

I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


Pig.  3 — The  bound  and  expected  distribution  of  the  normalized  correlation  of  independent 
random  sequences.  (The  expected  crosscorrelation  magnitudes  of  Gold-Pair  codes  are  shown 
for  comparison.) 


between  any  two  members  of  a Gold  family  of  codes).  The  upper  2 o values  are  equal  to 
the  mean  value  plus  twice  the  standard  deviation,  or  the  square  root  of  the  variance.  In 
the  case  of  random  sequences  with  Gaussian  distributions,  the  upper  2o  values  are  ex- 
ceeded with  a probability  of  only  0.025. 

Table  4,  Fig.  2,  and  Fig.  3 permit  a comparison  between  the  expected  properties  of 
the  magnitude  of  crosscorrelation  between  PN  Gold  pairs  and  pairs  of  purely  random 
sequences  in  relation  to  the  bound  determined  by  Gold  given  in  (58).  Several  maximal- 
length  sequences  of  period  L = 2N  - 1 are  examined,  where  N is  the  number  of  stages  in 
each  of  the  two  binary  shift-register  code  generators.  The  autocorrelation  of  independent 
random  codes  is  the  same  as  their  crosscorrelation.  Although  the  off-peak  autocorrelation 
magnitudes  of  Gold-family  codes  are  limited  by  the  Gold  bound  just  as  the  crosscorrela- 
tion magnitudes  are,  the  off-peak  autocorrelation  magnitudes  of  all  maximal-length  se- 
quences of  the  same  length  are  considerably  lower.  The  0-dB  bound  for  the  magnitude 
of  the  autocorrelation  or  crosscorrelation  of  random  codes  occurs  whenever  corresponding 
bits  in  the  two  sequences  happen  to  either  all  agree  or  all  disagree. 
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The  values  of  mean  magnitude  of  correlation  between  Gold  pairs  for  3 < N < 24 
are  from  about  0 to  4 dB  lower  than  the  corresponding  mean  values  for  independent 
random  sequences  (where  the  mean  is  taken  over  all  of  the  various  possible  code  phase 
conditions  to  obtain  the  expected  value  for  code  lengths  corresponding  to  those  produced 
by  N-stage  shift-register  code  generators).  The  expected  values  also  tend  to  be  as  much 
as  6 to  12  dB  below  the  Gold  bound.  The  superior  performance  of  the  Gold  codes  is 
due  to  their  especially  high  degree  of  mutual  orthogonality.  Adding  one  more  stage  to  a 
pair  of  odd-order  Gold  code  generators  allows  the  production  of  even-order  sequences 
with  a significantly  lower  mean  magnitude  of  correlation  but  with  only  a slightly  lower 
variance  and  bound,  whereas  adding  one  more  stage  to  a pair  of  even-order  Gold  code 
generators  provides  odd-order  sequences  with  a much  lower  magnitude  of  correlation 
bound  and  variance  but  only  a slight  improvement  or  reduction  in  the  mean  value.  An 
indication  of  the  standard  deviation  for  Gold  pairs  is  provided  by  the  lo  values,  which 
lie  above  the  mean  values  by  the  amount  of  the  standard  deviation. 

The  even-order  Gold  bound  is  comparable  with  the  upper  2o  values  for  independent 
random  codes;  hence  purely  random  sequences  can  be  expected  to  have  lower  correlation 
magnitudes  than  the  maximum  for  even-order  Gold  codes  of  the  same  length  for  about 
97.5  percent  of  the  cases.  However  odd-order  Gold  codes  have  a sharp  bound  (meaning 
that  they  can  at  times  achieve  the  bound)  about  2 dB  below  the  upper  2o  values  for 
random  sequences;  thus  independent  random  sequences  can  be  expected  to  have  larger 
magnitudes  of  crosscorrelation  than  the  maximum  for  odd-order  Gold  sequences  for  more 
than  2.5  percent  of  the  cases. 

The  primary  advantage  of  PN  Gold  codes  in  comparison  with  the  independent  ran- 
dom sequences,  in  addition  to  the  benefit  of  a lower  expected  value,  is  that  the  magnitude 
of  correlation  is  bounded  at  a significantly  lower  value,  which  decreases  even  more  as 
sequences  are  made  longer. 


SUMMARY 

The  magnitude  of  correlation  of  purely  random  sequences  has  been  evaluated  to 
provide  a fundamental  point  of  reference  needed  for  characterizing  the  relative  perform- 
ance properties  of  pseudonoise  sequences.  Expressions  for  the  mean  or  expected  value 
and  the  variance  for  the  magnitude  of  the  full-period  correlation  between  independent 
random  sequences  and  also  PN  Gold  pairs  have  been  determined  for  code  lengths  up 
through  those  corresponding  to  the  period  of  maximal-length  sequences  from  23-stage 
binary  shift-register  code  generators. 

The  comparison  of  the  crosscorrelation  characteristics  of  independent  random  se- 
quences and  PN  Gold  pairs  reveals  certain  advantages  of  Gold  pairs  due  to  their  inherent 
mutual  orthogonality. 

The  importance  of  choosing  longer  sequences  whenever  possible  is  clearly  indicated 
in  all  cases;  for  not  only  is  the  expected  magnitude  of  autocorrelation  and  crosscorrela- 
tion decreased  thereby,  but  the  variances  are  also  lowered.  Longer  sequences  thus  pro- 
vide better  correlation  performance,  whether  they  be  purely  random  or  pseudonoise  Gold 
codes,  achieved  by  the  use  of  more  stages  in  the  binary  shift-register  code  generators. 
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The  analysis  of  the  magnitude  of  correlation  for  random  sequences  provided  here  is 
useful  for  evaluating  the  expected  performance  of  signal  correlation  detectors  which  func- 
tion on  the  basis  of  the  magnitude  of  correlation  with  a reference  sequence  while  they 
are  operating  on  noise-only  input. 
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